Strong terahertz emission from superlattices via Zener tunneling 
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Abstract 

We develop a comprehensive, elegant theory to explain terahertz (THz) emission from a superlat- 
tice over a wide range of applied electric field, which shows excellent agreement between theory and 
experiment for a GaAsZAlo.3Gao.7As superlattice. Specifically we show that increasing electric field 
increases THz emission for low fields, then reduces emission for medium fields due to field-induced 
wave function localization, and then increases emission in the high field due to derealization and 
Zener tunneling between minibands. Our theory shows that Zener tunneling resonances yield high 
THz emission intensities and points to superlattice design improvements. 

PACS numbers: 73.21.Cd, 73.40.-c, 78.47.+p 
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The immense value of coherent terahertz (THz) electromagnetic radiation in many fields, 
including medical imaging, molecule recognition, sub-millimeter astronomy, remote sensing, 
and condensed matter physics^, has provided the impetus for vigorous research efforts into 
creating new sources of reliable and affordable THz radiation^ including exploiting Bloch 
oscillations in a superlattice (SL)4 The SL is especially promising as a THz source because 
it is effectively an artificially engineered semiconductor with tunable parameters and is 
thus adaptable. Vigorous experimental research^ 1 ^ has been accompanied by theoretical 
advances&^ii, but a proper theory that explains the physics and the THz emission properties 
for a SL THz source over a wide range of electric field has been elusive until now. 

We present a quantitative theory that explains THz emission from a SL over a wide range 
of applied electric field, and our theory also explains qualitatively the rise and fall of emitted 
THz field intensity in the low field (LF) regime, the rise of emission in the medium field 
(MF) to high field (HF) regime, anomalous resonances observed in THz emission from a 
GaAsZAlo.3Gao.7As SL£, and provides insight into the limitations in the ultrahigh field (UF) 
regime. Our theory elucidates the subtle physics of electromagnetic radiation emission in 
a SL subjected to high electric fields and, moreover, shows that Zener tunneling between 
minibands in the conduction band can be valuable especially for superior design of THz 
sources. 

Usually theories consider SL behavior in the regime of weak applied electric field where 
Zener tunneling (typically deleterious, such as electrical breakdown in semiconductors due 
to tunneling to higher bands^), can be ignored^. However recently we have seen an indi- 
cation that Zener tunneling between minibands of the conduction bancr^ coincide with the 
anomalous THz resonances observed in a GaAs/Al .3Gao.7As SlA Inspired by this promising 
connection between theoretical enhanced Zener tunneling rates and observed THz resonances 
in a SL, we have developed a theory for THz emission from a SL subjected to a range of ap- 
plied electric field; our theory explains experimental observations and provides a promising 
foundation from which to design SLs as versatile THz emitters. 

We analyze THz emission intensity / from a SL in an electric field F by subdividing 
the behavior into four regimes of applied electric field: (i) the LF regime corresponds to / 
increasing for increasing applied electric field F; (ii) the MF regime for which / decreases 
with increasing F; (iii) the HF regime for which I again increases with increasing F and also 
demonstrates resonances; and (iv) the UF regime (with I decreasing with increasing F) . Pre- 
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viously regimes (i) and (ii) have been separately explained by considering a single miniband 
within the conduction band using a Drude model for regime (i) and a quantum description 
for regime (ii)^; later the Kane model involving coupled Schrodinger equations^ was used 
to describe the cross-over between regimes (i) and (ii)^. More recently the locations of 
resonances in regime (iii) were identified with enhanced Zener tunneling rates between two 
minibands within the conduction bancU^ but the characteristics of / vs F were not pre- 
sented. Here we present a full and complete theory for regimes (i-iii) that includes all the 
successes of these prior models and, furthermore, easily describes all the cross-overs between 
applied electric field regimes, shows excellent quantitative agreement especially at the reso- 
nances. Our perturbative two-miniband theory provides a clear intuitive understanding of 
THz emission from a SL in terms of localization of wave functions and properties of both 
minibands. 

In our model the SL growth direction is along the x-axis with F in the positive x direction. 
For a and b the widths of the well and barrier, respectively, m the total number of unit cells 
(one well plus one barrier), d = a + b the unit cell width (and the potential V periodic 
over width d), and h the reduced Planck's constant, the Hamiltonian for an electron — e and 
effective mass \i is 

h 2 d 2 

H = - — — + V(x) + eFx,V(x) = V(x + ld). (1) 
2/i ax z 

The spectrum is subdivided into minibands indexed by n, and we are interested in the 
case of two minibands so n e {1,2}. The eigenfunctions for F = can be expressed as 
Bloch functions {4> n k{x)}, with corresponding eigenvalues {E n k} and Bloch number k in 
the Brillouin Zone (BZ). The Bloch functions and the eigenvalues in the SLs are obtained 
by using the Kronig-Penney modet^ with a step size of d/200 and the BZ is sampled by 
2000 points. In order to appreciate the role of inter-miniband tunneling, we study the 
system with Zener tunneling between minibands treated by perturbation theory; exact and 
numerical methods for solving Hamiltonian H (JT]) tend to obscure the important role played 
by tunneling. Furthermore we will see that the ultrahigh electric field regime, which is 
characterized by decreasing J, is where second-order perturbation effects become important. 

The Bloch functions can be used as a basis set for F ^ as well. In this case, a 
wave function confined to band n is $ n ;(x) = XlfceBZ a nik4>nk{x) , with the coefficients {a n ^} 
satisfying a well known relation^^. Inter-miniband tunneling is calculated via dual degen- 
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erate perturbation theory. The zeroth-order wave function V'i°2/'( x ) = Cu^u(x)+C2i'^2i'(x) 
straddles both minibands, with V satisfying the condition that the energy difference | e 2 r — £u\ 
is minimized in V G {1, . . . , m}, and 
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with £y the eigenvalue for first-order perturbation of Eq. ([1]) and H' u 2l , = 
eFx | §2i'{x)) the perturbation matrix element. The first-order wave function 

is 

/W ( \ _ ^i«",2r-«+i"( :r )l e -^ x I ^1^,2/' C^-) ) 

VU,2l'\ X ) - 2^ (0) _ (0) 

£ ll £ ll" 

X ' l l J ll'',2l'~l+l"( X )i (3) 

with I" G {1, • • • and the zeroth order eigenvalue. In the case of weak tunneling 
(Cu — > 1, — > 0), degenerate perturbation theory merges with non-degenerate pertur- 
bation theory. For first-order perturbation theory, the total wave function is ipu,2i'{x) = 
'tPu^i' ( x ) + V'iI \vi x )- We also perform the calculation to second order in order to separate 
first-order effects, which yield desirable resonances, and second-order effects, which we show 
limit performance in the ultrahigh-field regime. Whether the wave function ipu t 2i'(x) is 
calculated to first- or second-order will be clear from the context. 

We plot Re[ipn^i'(x)} in Fig. [T] for several values of F to depict characteristics of this 
wave function such as localization and overlap with the nearest levels in each of the four F 
regimes (LF, MF, HF, UF). We see that the wave function can be localized or delocalized 
and may be confined to one miniband or overlap two minibands, with these characteristics 
important for understanding THz emission. 

THz emission intensity depends on the wave function ^11,2/' (x) and on the applied electric 
field. The relationship between THz emission frequency ujb and field strength F is given 
by hjjB = eFd, which is the energy level splitting in the miniband. To calculate the total 
radiative intensity, we adapt the relation^ 



m „ 

fJ2\i-i"\ 
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dx ^11 21' (x)xip U ",2l'~l+l" (x) 
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'SL 

to include broadening due to scattering; here / G {1, . . . , m}, V minimizes \e 2 °J — e$\, and 
the wave functions implicitly depend on F. For the experiment to be considered in the 



present paper, the energy half width of the pumping laser pulses is approximately 20 meV. 
Within this energy interval, the levels in the Wannier-Stark ladders (WSLs) are populated 
with equal probability^. Thus the summation runs over all pairs of WSL eigenfunctions. 

It has been demonstrated in Ref.— ^ that the acoustical phonon scattering plays a decisive 
role in the THz radiation process. Let \n, p) be the eigenstate of the nth WSL level, where 
p is the two-dimensional (2D) wave vector. As shown in Ref.— , in the steady state the 
electron population on each WSL level is the same, and there will be no net stimulated 
emission from the SLs due to the vertical transition from the initial state \n + l,p) to the 
final state \n, p) is perfectly canceled by the absorption transition from state \n, p) to state 
\n + l,p) that occurs at the same frequency and intensity without the effect of electron- 
phonon scattering^. However, in the transition processes involve the emission or absorption 
of an acoustical phonon, the in-plane momentum of the initial state will differ from that of 
the final state. Hence, the emitted THz radiation will not be completely reabsorbed by the 
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FIG. 1: (Color on-line) Plot of Re [ipnpi' {%)] (blue line), which is the wave function at energy I, I' for 
applied electric field values (a) F = 14kV/cm (LF), (b) F = 26kV/cm (MF), (c) F = 45kV/cm (HF 
away from THz resonance), (d) F = 90kV/cm (HF near THz resonance), and (e) F = 120kV/cm 
(UF). The horizontal lines are the equally spaced energy levels in each of the two minibands (right 
for the 1st miniband and left for the 2nd miniband). 
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FIG. 2: (Color on-line) First-order (solid line) and second-order (dashed line) perturbation theory 
results vs experimental data [8] (black crosses) for structure 1. 

system, and the system can have a net gain for the photon energy fiw < Hub- Because the 
energy of an acoustical phonon is much less than the energy of THz photon, we still treat 
the hwB as the THz radiation photon energy in this work-. 

Our theory is compared with recent experimental results on two GaAs/Alo.3Gao.7As SL 
structures^: structure 1, a = 6.4nm and b = 0.56nm with m = 73 (first miniband: 18 to 114 
meV; second miniband: 150 to 445 meV) vs structure 2, a = 8.2nm and b = 0.8nm with 
m = 55 (first miniband: 19 to 69 meV; second miniband: 107 to 270 meV). We show our 
calculated I{F) characteristic for both first- and second-order perturbation theory in Figs. [2] 
and [3] alongside corresponding experimental results for both structures. 

Two fitting parameters are used in our calculations to match theoretical and experimental 
curves: If lot (F) = ^qI^IF) + <^i, with If lot the plotted / for structure £, i| alc the calcu- 
lated /, and the fitting parameters f° r rescaling and for background. The choice of 
both coefficients is jointly determined by least squares fitting for both sets of experimental 
data, with excellent agreement between second-order perturbation theory and experiment. 
First-order perturbation results agree for all electric field regimes except UF, thus revealing 
the important role of second-order effects in limiting achievable THz intensity. 

Theoretical and experimental radiation intensities I{F) are depicted in Figs. [2] and [3] for 
structures 1 and 2, with the theoretical curve determined by first- and second-order pertur- 
bation based on Eq. (jlj). As expression (jlj) ignores scattering, we incorporate the broadening 
effects by convolving 1(F) with a Lorentzian of half-width T, which is obtained based on 
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FIG. 3: (Color on-line) First-order (solid line) and second-order (dashed line) perturbation theory 
results vs experimental data [8] (black crosses) for structure 2. 

TABLE I: Fitting parameters for structure 1 and 2 



Parameters 




ft 


Structure 1 (lst-order) 


3.667 


0.316 


Structure 1 (2nd-order) 


3.739 


0.332 


Structure 2 (lst-order) 


2.598 


0.164 


Structure 2 (2nd-order) 


2.771 


0.134 



the calculation of longitudinal optical phonon scattering^^i^, aluminum atom alloy scat- 
tering^ 1 ^ and interface roughness scattering^ 1 ^ by using Fermi's golden rule. Theory and 
experiment agree well for all field regimes except UF, and first-order perturbation results 
also agree for all applied field regimes except the UF. However note that the theoretical 
curve for structure 2 diverges somewhat from the UF intensity peaks of the experimental 
curve, which we believe is due to these two peaks arising from Zener tunneling between the 
second and third minibands; furthermore this anomaly is compounded by our least square 
fitting parameters So, and ft truncated for F cut = 100 and F cut = 120 kV/cm for structure 1 
and 2, respectively. The fitting parameters for the first-order and second-order perturbation 
calculation in structure 1 and 2 are given in table E 

The curves clearly exhibit four regimes of behavior according to the strength of the applied 
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electric field. We see that 1(F) rises monotonically for LF, then decreases monotonically 
in MF, and HF is characterized by a trend of / increasing with rising F but punctuated 
by resonances. Finally there is a region in the UF case with the experimental result for I 
decreasing with increasing F, which agrees with the second-order perturbation theory result, 
but disagreeing with first-order perturbation result. Our theory explains all these features 
as exemplified by Figs. HJ In Figs. l(a,b), corresponding to the low- and medium-field 
regimes, respectively, the slopes of the minibands are small, and Zener tunneling between 
minibands can thus be ignored. In these two cases the wave functions are confined to 
the lowest miniband, but in (a) the function progressively spreads out with increasing F 
whereas it progressively localizes with increasing F (due to competition between localized 
wave functions) in (b)^ 1 ^. 

The HF regime in Figs. QJc-d) exhibits clear effects due to inter-miniband tunneling. 
Fig. 1(c) shows steeper minibands than in[Tfa-b), which causes the wave function to overlap 
both minibands. Analogous to (a), the wave function progressively spreads as F increases, 
yielding increasing THz emission with increasing F . However this increasing / with increas- 
ing F is punctuated with resonances due to Zener tunneling (this connection was speculated 
when coincidences between values of F values for enhanced Zener tunneling and experimen- 
tally observed resonances were demonstrated^). 

Here we see exactly what is happening at the resonance: the wave function is localized 
in both minibands but with significant overlap, which is helped by alignment of the WSLs 
in each miniband, and hence are strongly coupled leading to enhanced THz emission. In- 
creasing F partially leads to further localization, as in the medium-field regime, but causes 
the WSLs to become misaligned, thereby decreasing the Zener tunneling rate and causing a 
decrease of / with increasing F. In this way, we can understand the existence of these THz 
resonances, and, moreover, design SLs to exploit these emission peaks; moreover our theory 
makes it evident that, when second-order perturbation events become non-negligible, wave 
functions become progressively delocalized thereby making the SL an increasingly poor THz 
emitter. 

In summary we have developed a clear model of electromagnetic emission from a SL, 
based on single miniband dynamics in the low- and medium-field regimes, and two-miniband 
dynamics in the high- and ultrahigh-field regimes. Our results contain previous theories for 
low- and medium-field dynamics as special cases and explain the coincidence of enhanced 
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Zener tunneling rates and experimentally-observed THz emission peaks. We have excellent 
quantitative agreement with experiment, subject to choosing two fitting parameters by using 
least squares fitting method, and we have a clear explanation of when and why THz emission 
becomes poor in the ultrahigh applied electric field regime. In addition to the value of our 
model as a design tool for SL THz emitters, our perturbation theory approach is remarkably 
successful in explaining electromagnetic radiation features in recent SL experiments due to 
confinement or coupling between minibands within the conduction band. Our approach also 
underscores the value of a perturbation theory approach, as opposed to exact or numerical 
means, to understand the underlying physics. 
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